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1. Introduction 

The study of heat kernels is an important problem in the theory of parabolic 
PDE’s. The properties of heat kernels give an efficient tool to answer to 
some central questions both in analysis and probability theory. One of the 
main questions is to know whether a heat kernel admits Gaussian bounds. 
An upper Gaussian bound is for instance an useful tool for getting 
estimates, the analyticity of the corresponding semigroups in for any finite 
p > 1 or bounded functional calculus, whereas one can get a strong maximum 
principle or a Harnack inequality from a lower Gaussian bound. We refer to 
the textbooks 0: E] and m and references therein for more details on the 
subject. 

In the preceding work [3] , starting from the classical parametrix method, 
we constructed the Neumann heat kernel of a general parabolic operator as 
a perturbation of the fundamental solution of the same operator by a single¬ 
layer potential. From this construction, the two-sided Gaussian bounds for 
the fundamental solution and taking into account the smoothing effect in 
time of the single-layer potential, we succeeded in proving a lower Gaussian 
bound for the Neumann Green function. We adapt in the present note this 
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method to establish a lower Gaussian bound for the Neumann heat kernel of 
Laplace-Beltrami operator. 

In this text A4 = (A4 ,g) is a n-dimensional compact connected Rie- 
mannian manifold without boundary and is a domain in j\4 so that its 
boundary E is an (n — l)-dimensional Riemannian submanifold of M when 
it is equipped with the metric induced by g. 

The Riemannian measure on M is denoted by dV while the density 
measure on E is denoted by dA. The geodesic ball of center x £ M and 
radius r > 0 is denoted by B{x,r). 


2. Neumann heat kernel 

Let d be the Riemannian distance function and 

(S'{x,yA) = (47rt)“”'^^e . 

In general S is not a heat kernel on M.. However the parametrix method 
by S. Minakshisundaram and A.Pleijel shows that any compact Riemannian 
manifold has an almost euclidien heat kernel (e.g. BE])- In particular, the 
heat kernel p oi M., satisfies 

P{x,y,t) ^ S‘{x,y,t), 

locally uniformly in (x, y) as t 4- 0- We have a similar statement with the first 
order derivatives of p and S. 

This estimate is not true in general for distant x and y. Indeed the 
counter example in [M Example 3.1, page 23] shows that if ^ — S^, the 
2-dimensional sphere equipped with the round metric, x is the north pole 
and y is the south pole, then 

P{x,y,t) ^ 

for some constant c > 0. 

From 0 Theorem 5.5.II and Theorem 5.6.1, page 173], any complete 
Riemannian manifold with non negative Ricci curvature satisfies the following 
two-sided Gaussian bounds. 

<S’{x,y,t) <pix,y,t) < £ M, t > 0. 

V[B{x,yyt)) 

Here c > 0 and k > 0 are some constants. 

Let A = Ag be the Laplace-Beltrami operator associated to the metric 
g and denote by v the outward normal vector field to E. Following the idea 
in [3], we construct the Green function of the Neumann problem 

f (dt — A)u = 0 in H X (0,-l-oo), 

I If = 0’ on E X (0,-foo), 


( 2 . 1 ) 
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as a perturbation of the heat kernel p by a single-layer potential. As a first 
step, we seek the solution, in x (0, -koo)) fl x [0, -koo)), of the 

following IBVP 

{dt — A)m = 0 in n X (0, -koo), 

fH = 0, onSx(0,+oo), (2.2) 

u(-,0) = V^GCo°°(r!) 

of the form 

u{x,t) = / / p{x,y,s)ip{y,t-s)dA{y)ds+ / p{x,y,t)ip{y)dV{y), 

Jo Jy. Jn 

{x, t) G n X (0, -koo). 

We obtain from the jump relation in [2 Theorem 2, page 161] that (p must 
be the solution of the following integral equation 

p{x,t) =-‘J [ f ■^ix,y,s)p(y,t-s)dA{y)ds 


0 JY 


dvx 


- 2 


f ^ 

In dvx 


[x, y, t)ip{y)dV(y), (x, t) G S x (0, -koo). 


This integral equation is solved by successive approximations. We get 

(p{x,t) = Poix,t)+ / r{x,y,s)ipo{y,t-s)dA{y)ds, (x, t) G S x (0,-koo). 
Jo Jt. 


Here 


ipo{x,t) =-2 [ -^ix,y,t)'tpiy)dV{y), (x, t) G S x (0,-koo), 

Jn 

r{x,y,t) = '^r^{x,y,t), ix,t) G E x (0,-koo), y G fl, 

j>i 


with 


dp 


r ix,y,t) = -2 — {x,y,t), {x,t) e S x (0,-koo), y G fl, 

Ol^x 

r^+^{x,y,t) =-2 / r^{x,z,t-sy{z,y,s)dA{z)ds, j > 1, 

Jo Jy 


(x, 1) G E X (0, -koo), y G fl. 


In other words 
, dp 
dvx 


r{x,y,t) =-‘2'-^ix,y,t)-2'^ f f {x,z, s)-^{z,y,t - s)dA{z)ds 

0 J ^ X 


(x, 1) G E X (0, -koo), y G fl. 

The following two inequalities will be useful in the sequel. They are 
taken from [2. 

For any /r > 0, there is a constant Co such that 

dp 


dv. 


-{x,y,t) 


<Col-'^d-”+"^(x,y), 


(2.3) 
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for any x G Y,, y G and t G (0,T]. 

There exists a constant Ci such that, for any a, ^ G (0, n — 1), 

if Of +/3 > n — 1, 

1 ifa + /3<n — 1. 

(2.4) 

for all x,y G X y. 

With the help of inequality (I2.d|l . we prove similarly to O (3.6)] 


d °‘{x,z)d ^{z,y)dA{z) < Cl 


(p{x,t) = / r{x,y,t)ipiy)dV{y). 
Jn 


Then 


where 


u{x,t) = / qix,y,t)'ipiy)dV{y), 

Jn 

qix,y,t) = Pix,y,t) + / / pix,z,s)r{z,y,t-s)dA{z)ds. 

Jo Jn 

We call this function the Neumann heat kernel for the problem (EU. We leave 
to the reader to verify that q satisfies the following reproducing property 

q{x,y,t)= qix,z,t-s)q{z,y,s)dv{z), x,yG^,0<s<t, 

Jn 

and 

[ q{x,y,t)dV{y) = 1, xGn,t>0. 

Jn 


3. Gaussian lower bound 

We set v{x,r) = V{B{x,r) fl fl), x G ft, and we consider the following three 
assumptions. 

(VLB) (volume lower bound) There exist two constants C and tq so 

that 

v{x,r) > Cr"^, X G fl, 0 < r < tq. 

{DP) (doubling property) There exist two constants ri > 0 and C > 0 
so that 

( s \ ^ 

-j v{x,r), 

for all 0 < r < s < ri and x G fl. 

(CC) (chain condition) There exists a constant C > 0 such that for any 
x,y G fl and k G N, we find a sequence of points {xi)o<i<k so that xq = x, 
Xk = y and 

d{xi,Xi+i) < C ^ 0<i<k — l. 

k 

In the flat case, (VLC) and (DP) are true for any Lipschitz domain 
and (CC) holds for instance for a convex domain. 
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Let be the unit sphere of equipped with the round metric, that 
is the metric induced by the Euclidean metric on Then any segment of 
length strictly less than tt possesses the the three conditions {VLB), {DP) 
and {CC). Let be the unit sphere of equipped with the round metric. 
Then it is not hard to show that the sub-domain of given by 0 > <5 > 0 
satisfies also the three conditions {VLB), {DP) and {CC). 

Conditions {DP) and {CC) are usual (see for instance [HI Theorem 
7.29, page 248]), while condition {VLB) guarantees that the near diagonal 
lower bound (13.101) above holds. 

We aim to sketch the proof of the following theorem 

Theorem 3.1. Fix T > 0 and assume that 11 satisfies {VLB), {DP) and 
{CC). Then 

q{x, y, t) > ^ ^ e ^, x,y en, 0 <t<T. (3.5) 

v{x, y/t) 

Sketeh of the proof. Let 1/2 < p < n/2. In light of (I2.3|l and (12.4L reasoning 
as in the proof of [Sj Lemma 3.1], we obtain 

\r{x,y,t)\<Ct-^d-^+^^{x,y), 

for any x € T,, y € fl, x ^ y, t € (0,T]. 

Let 

qo{x,y,t) = / / p{x,z,s)r{z,y,t-s)dA{z)ds, x,yGLl,t>0, 

Jo je 

and 0 < a < 1/2. We proceed as in the beginning of the proof of [5] Theorem 
3.1] to get 

\qo{x,y,t)\ < x,yGn,t>0. (3.6) 

Let e := inj(Al)/4, where inj(A4) is the injectivity radins oi A4. It 
follows from [31 formula (45), page 154] that there exists y > 0 so that 

p{x,y,t)> <g{x,y,t), 0 <t<r], x,y€ J^, d{x,y) <e. (3.7) 

Hence, 

p{x, y, t) > Ct~'^l'^, 0 < t < mi{r], e^), x,y € J^, d{x, y) < y/t. (3.8) 
Now a combination of (13.611 and (13.811 leads 
q{x,y,t) > p{x,y,t) - \qo{x,y,y)\ > Ct~'^l‘^{\ - cC), 

0 <t < mf{ri,e^), x,y € fl, d{x,y) < y/t. 
In consequence, there is <5 > 0 such that 

q{x,y,t) > Ct~^^^, if 0 < t < <5, x,yGLl, d{x,y)<y/i. (3.9) 

In light of the volume lower bound {VLB), this estimate entails 
C ~ 

q{x,y,t)>— -if 0 < < < x,y G Q d{x,y) < y/i, (3.10) 

v{x, y/t) 


for some constant S. 
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We can now mimic the proof of El Theorem 7.29, page 248]. We get 
from (|3.10l) the following Gaussian lower bound 

Q d^(x,y) 

q{x,y,t) > — - x,y en, 0 <t<S. 

v{x,Vt) 

We finally use the argument as in [21 Theorem 3.1] to pass from 0 < t < S to 
0 <t<T. □ 

It is worthwhile mentioning that one can establish a lower Gaussian 
bound by considering Af = Q itself as a compact Riemannian manifold with 
boundary. The structure of Riemannian manifold is the one inherited from 
M. Obviously, (1231) entails 

q{x,y,t) > , 0 < t < S, x,y€j\f, djq-{x,y) < Vi. (3-11) 

Assume that the Ricci curvature of Af is such that Ric > {n — l)Kg for some 
K gM.. Then Af satisfies the a doubling property (DP) when V is substituted 
by V/v", the volume measure over Af. This fact is an immediate consequence 
of [71 formula in the bottom of page 7]. Let vjq-{x,r) = V_\f{B{x,r)), where 
B{x, r) is the geodesic ball in Af of center x G Af and radius r. In that case we 
can proceed as in the proof of Theorem 13.11 to derive the following estimate. 

g(a;, y, t) > — ^ e , x,y Gn, 0 <t<T. 

v^ryx, y/t) 

This estimate should be compared to the one obtained by Li and Yau in [21 
Theorem 4.2, page 184]. Specifically, they established a lower Gaussian bound 
for the heat kernel of a compact Riemannian manifold with convex boundary 
and having non negative Ricci curvature. 

When AA is any complete Riemannian manifold with finite diameter 
and having volume doubling property, and is Lipschitz domain in AA with 
volume doubling property, the Neumann heat kernel of 12, denoted here by 
h, satisfies the following upper Gaussian bound. 

C d^(x,y) 

Hx,y,t)<— -r=—- —e 8* , x,y Gn,t>0, 

v{x, Vt)v{y, y/tj 

where C > 0 is some constant. 

This estimate was recently established by the authors and E. M. Ouhabaz 

i- 

4. Comments on geometric assumptions 

Chain condition: A subset C of A4 is called strongly convex if for any x,y G C, 
there exists a unique minimal geodesic 7 : [0,1] —> A4 joining x to y, so that 
7([0,1]) C C. According to a theorem due to Whitehead (see for instance [21 
pages 161 and 162]), there exists a positive continuous function e : A4 —>■ 
(0, 00], the convexity radius, such that any open ball B{x,r) C B{x,e{x)) is 
strongly convex. It is straightforward to check that if 12 is strongly convex 
then it has the chain condition. 
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Volume lower bound: Let T^Ai be the tangent space at a: S C 

TxM the unit tangent sphere and SAi the unit tangent bundle. Let be 
the geodesic flow with phase space SAi. That is, for any t>0, 

<i>t : SAi ^ SA4 : {x,0 G SAi ^ • 

Here 'yx,^ : [0, oo) —>■ is the unit speed geodesic starting at x with tangent 
unit vector ^ and 7a;,^(t) is the unit tangent vector to at 7a;,{(t) in the 
forward t direction. 

If {x, S SAi, we denote by r(x, the distance from x to the cutlocus 
in the direction of ^: 

r(x, = inf{t > 0; d(x, $t(x, ^)) < t}. 

We fix i5 S (0,1] and r > 0. Following [T3], a (<5, r)-cone at a; G AI is the 
set of the form 

‘^{x,ujx,r) = {y = 7x,i(s); ^ € lux, 0 < s < r}, 

where uix is a subset of Sa; so that r < r{x,^) for all ^ G Wx and |a;a;| > 5 
(here \liJx\ is the volume of uix with respect to the normalized measure on the 
sphere Sa;). 

A domain D which contains an (6, r)-cone at x for any x G D is said to 
satisfy the interior (<5, r)-cone condition. 

We observe that if C is a closed strongly convex subset of AI, then 
id = Ai\C has the (1/2, r)-cone condition, for some r (this fact follows from 
the same argument to that in [T31 Example 8.1, page 370]). 

Let 

f (=^)”’‘ if-x). 

Sk(?') = < r” ^ if K = 0, 

[ ifK<0. 

We make the assumption that the sectional curvature of Ai is is bounded 
above by a constant k, k G M, and fl satisfies the interior (6, r)-cone condition. 
Let J{x,^,t) be the density of the volume element in geodesic coordinates 
around x: 

dV{y) = J{x,^,t)ds,,dt, y = 7x,^{t), t < r{x,Cj. 

By an extension of Gunther’s comparison theorem (see for instance [S]), J 
satisfies the following uniform lower bound 

> SK(t). 

Consequently, shrinking rg if necessary, we have 

v{x,r) > V{^{x, 0 Jx,r)) > co'r"', x G 0 < r < ro, (4.12) 
which means that v satisfy the volume lower bound {VLB). 

Additionally, if Al satisfies the following volume growth condition 

H(a;,r) < cir", 0 < r < tq, (4.13) 
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for some constants ci and ri then v has the doubling property {DP). 

As a consequence of Theorem 13.11 we have 

Corollary 4.1. Assume that the sectional curvature of M is bounded from 
above, the volume growth condition (14.131) is fulfilled and Vl is strongly convex 
and satisfies the interior {5, r)-cone condition. Then 

q{x, y, t) > cS'{x, y, ct), x,y € D, 0 < t < T. (4.14) 
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